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■ In this paper we study the gluon distribution functions in nucleons and pions 

at a low resolution scale. This is an important issue since parton densities at 
low have always been taken as an external input which is adjusted through 
DGLAP evolution to fit the experimental data at higher scales. Here, in the 
framework of a model recently developed, it is shown that the hypothetical 
cloud of neutral pions surrounding nucleons and pions appears to be responsible 
for the characteristic valence-like gluon distributions needed at the inital low 
scale. As an additional result, we get the remarkable prediction that neutral 
and charged pions have different intrinsic sea fiavor contents. 
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1 Motivation 



Parton distribution functions (pdf) play an important role in particle physics as they 
describe the internal structure of hadrons in the framework of Quantum Chromody- 
namics. Pdf are also basic ingredients to calculate cross sections in hadron-hadron 
and lepton-hadron interactions. 

Parton distribution functions have unique characteristics depending on each ha- 
dron, which reflect the internal dynamics of the bound state. The pdf, are 
interpreted as the probability of finding a parton i (quark or gluon) with a fraction 
X of the hadron momentum when probed by the momentum transfer Q^. The 
dependence of pdf is succesfuUy described by the DGLAP |I| evolution equations 
within perturbative QCD. However, in order to have adequate fits to Deep Inelastic 
Scattering (DIS) data, initial valence-like distributions at a low resolution scale must 
be considered even for sea quarks and gluons. This fact has been repeatedly noted by 
several authors (see e.g. Refs. p|, 0, H) but the origin of such initial distributions 
still remains unclear. 

The origin of primordial low scale pdf should be traced back to the internal dy- 
namics of the hadron bound state. Thus they must be related to the confining phase 
of QCD, which is ultimately the origin of hadrons as bound states of quarks and glu- 
ons. In this sense, the sea quark and gluon distributions at a low resolution scale can 
be related to the idea of the intrinsic sea of quarks and gluons proposed by Brodsky 
et al. at the beginning of the 80's |^. In their own words, "intrinsic quarks and 
gluons exist over a time scale independent of any probe momentum, and are associ- 
ated with the bound-state hadron dynamics". In contrast, the extrinsic sea of gluons 
and quarks have a purely perturbative origin, and their distributions show the char- 
acteristics of bremsstrahlung and pair production processes leading to the standard 
DGLAP perturbative QCD evolution. 

The above considerations led us to draw a low picture of hadrons in terms 
of effective quark degrees of freedom interacting among them through intrinsic sea 
quarks and gluons . The building blocks are the so-called valons which are valence 
quarks appropriatedly dressed by their extrinsic sea Within this picture one can 
represent the hadron wave function as a superposition of hadron-like Fock states which 
we construct by means of a well-known recombination mechanism For example, 
the proton wave function can be written as 

\P) = flO bo) + ai \po g) (1) 

i 

at some low Ql scale compatible with the valon picture of the proton. Here \po) is a 
pure three valon state 0, and the following terms are hadronic quantum fluctuations 
which emerge from the (non-perturbative) interaction among valons. The role of 
these hadronic fluctuations is to dynamically generate the intrinsic sea of quarks and 
gluons, which is the necessary binding agent in order to have the hadron state. In 
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this way, a consistent picture of the low (input) scale of hadrons emerges in which 
hadrons are formed by constituent quarks plus intrinsic sea quarks and gluons. 

It is worth to stress that in modern fits to DIS data, initial non-perturbative sea 
quark and gluon distributions are taken as an input which is adjusted by DGLAP 
evolution. In contrast, in our approach they are dynamically generated through the 
hadronic quantum fluctuations. Furthermore, the old fashioned fits to DIS data, in 
which all the hadron sea is perturbatively generated, are recovered by restricting the 
series of eq. (|I|) to the first term. 

This representation of the proton wave function led in a natural way to a (i/w 
asymmetry in the proton sea which closely describes [|T0| the most recent experimental 



data by the E866/NuSea Collaboration On the same footing, in Ref. [0] a s 



asymmetry in the nucleon sea was calculated, qualitatively agreeing with the results 
of the last global analysis of DIS data fT^ . 

In the following section of this paper, we determine the low non-perturbative 
gluon distributions in nucleons, using the model introduced in Ref. . In section ^ we 
explore the consequences of this picture on the non-perturbative structure of charged 
and neutral pions. Finally, section H is devoted to further discussion and conclusions. 



2 Intrinsic gluon distributions 

To start with, let us consider a baryon at some low Q\ scale. At this scale the baryon 
ground state is formed only by three valons 0. Quantum fluctuations will generate 
the non-perturbative qq sea. Following Ref.0, the non-perturbative sea has a two- 
step origin in our model. In the first step a valon emits a gluon which subsequently 
decays into a quark- ant iquark pair. In the second, such quark and antiquark interact 
with the valons giving rise to a bound \MB) state. Non-perturbative quark and 
antiquark distributions are then associated to the in-hadron meson and baryon valon 
densities. 

The emission of a gluon out of a valon is a basic QCD process which can be 
adequately described in terms of the convolution of the valon distribution, v(^z\ with 
the Altarelli-Parisi Pgq{z) and Pqg{z) splitting functions |I|]. In this way, the quark 
and antiquark initial distributions are given by 

Notice that, as the valon distribution does not depend on Q^, the scale dependence 
in eq. only arises through the strong coupling constant Ugt- At this stage the 
scale is fixed to the valon scale, which is tipically about 1 GeV^. Indeed, in 
Ref. the valon scale was estimated to be of the order of Q\ ~ 0.64 GeV^ for 
nucleons. Consequently, pair creation can be safely evaluated in a perturbative way 
since (asj/27r)^ is still sufficiently small. 
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The perturbative v ^ g qq process is the source of both the extrinsic and 
intrinsic seas. The difference among them rest on the fact that an intrinsic qq pair 
interacts with the remaining valons while an extrinsic qq pair not. In this sense, the 
extrinsic sea, which is purely perturbative, forms the structure of valons 0. 

The second step involves the interaction of such qq pair with valons thus giving rise 
to the \MB) bound state. As they are in the realm of confinement, the interactions 
of the qq pair with valons must be evaluated by means of effective methods. Notice 
also that, for such interactions take place, the initial (perturbative) qq pair must be 
sufficiently longlived. Since the characteristic lifetime of such a perturbative qq pair 
scales as l/rriq, light and strange quarks should be largely available to interact with 
valons thus producing the \MB) hadronic quantum fluctuations. 

Then, assuming that the in-hadron meson and baryon formation arises from mech- 
anisms similar to those at work in the production of real hadrons, we can evaluate 
the in-hadron meson probability density by using the Das-Hwa recombination ap- 
proach 0. 

In the recombination model, the probability density for the production of a real 
meson as a function of its fractional momentum is given by the convolution of a 
two-quark distribution with a suitable recombination function. The two-quark distri- 
bution is given in terms of the single-quark distributions of the initial hadron which 
will be the valence quarks in the final meson. The recombination function is chosen 
in such a way that it favors the recombination of quarks with similar momentum 
fractions. 

Thus, in our model, the in-hadron meson distributions are given by 

/■I (^y rl 

Pm,bX^)= — —Ft{y,z)R{x,y,z), (3) 
Jo y Jo z 

where 

= ap5(l-^) (4) 
is the recombination function 0, and 

Fi{y,z)=Pyv{y)zq,{z){l-y-zr (5) 

is the valon-antiquark distribution |TD||. In eqs. (^-(^, x, y and z are the momentum 
fraction of the in-hadron meson, the valon and the antiquark respectively. The index 
i runs over different quark flavors, depending on the meson being formed. 

Due to momentum conservation, the in-hadron meson and baryon probability 
densities are not independent but correlated by 

Pm^bX^) = Pb^mX^ - , (6) 
with an additional correlation in velocity given by 

{xPm^bA^)) ^ {xPb^mA^)) /r,x 
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The above constraint, eq. (J^), which is needed in order to build a \MiBi) bound state, 
fixes the exponent a in eq. [p!0|] . 

The hadronic fiuctuations so far described can be interpreted as the origin of 
the intrinsic quark-antiquark sea. As a consequence, since the resulting q and q sea 
distributions belong to different hadronic states in the \MiBi) fluctuation, intrinsic 
quark and antiquark probability densities in baryons are unequal in a general way. 

At this point, a judicious analysis of what fluctuations should be included in an 
expansion like eq. (|lD must be made. For deflnitness, consider the proton wave func- 
tion. Taking into account mass values and quantum numbers, the main fluctuations 
of the proton should be the n), \7i~^ A^) and |7r~ virtual states, with proba- 

bilities |ci7r„P and |a^Ap respectively. Differences between the |7r+ A*^) and |7r^ 
probabilities are taken into account by Clebsh-Gordan coefficients which ensure the 
correct global isospin of the fiuctuation. Thus, we obtain IIottaP and ||a^Ap for 
Itt"*" A") and |7r^ A"^"*") respectively. On the other hand, the probability of the \7r~^ n) 
bound-state is ||a7rnP- The coefficients lajmp and IottaP are given by {NaP)nn and 
(A^a/3),rA respectively. Their numerical values result from comparison with experi- 
mental data. Fluctuations like \pN), |pA), etc. which could contribute, for instance, 
to the d — u asymmetry in the proton, are far off-shell Q and can be safely neglected. 

Remarkably, as shown in Ref. this scheme leads to a d/u asymmetry in the 



proton which closely describes the experimental data by the E866 Collaboration ||TT 



Including fluctuations to Kaon-Hyperon states, \K H), a s — s asymmetry in the pro- 
ton sea arises which qualitatively agrees with results from a recent global analysis 
of DIS data 0. 

Fluctuations of the proton to \tt^ p) and Ivr" A+) states do not contribute to the 
intrinsic quark and antiquark structure. The reason is that the formation of a 7r° 
in-proton state must be inhibited due to its neutral flavor structure uu — dd. This 
happens because Vgq objects can annihilate rapidly into a gluon while Vgq' cannot [q' 
is a quark of different flavor to the q one. See Fig. |^). Notice also that an unflavored 
object like avqq pair has itself the quantum numbers of a gluon. Thus, a hypothetical 
|7r°p) fluctuation does not contribute to the sum over \MiBi) in the RHS of eq. ([^) but 
to the second term, \po g), providing a source of valence- like gluons in the proton. The 
proton-like object accompanying the non-perturbative gluon in the \po g) fluctuation 
must have the same flavor structure than the Po in the \po) state. However, as far as 
the gluon is in a color octect state, the po must be colored. It is worth noting that, 
on general grounds, hadrons in a \MB) fluctuation must be colored. However, they 
can be identified with real hadrons regarding other quantum numbers like as fiavor, 
isospin, etc. 

The time scale over which both the \n~^n) and the \po g) fiuctuations exist should 
be approximately the same. In fact, the characteristic lifetime of a \M B) fiuctuation 
is proportional to 1/ AE, where AE is the energy difference between the \M B) and 

^Note that these are even more suppressed than strange, \K H), fluctuations. 
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the proton states in an infinite momentum frame. Tlius, for a generic \M B) state we 
have 

1 2P 



'^\MB) 



AE 



B _ 



(8) 



where P is the momentum of the proton in the infinite momentum frame, rup is the 
proton mass, and xm and xb are the momentum fractions carried by the meson and 
baryon in the fluctuation. m\j ^ = ^"ii b + is the transverse masses squared of 
virtual hadrons in the fluctuation. Given the smallness of the pion mass, we can 
assume that in-nucleon pions and non-perturbative gluons have similar transverse 
masses, then the characteristic lifetimes of the Ivr"*" n) and \po g) fluctuations must be 
approximately the same. 

In this approach, the shape of the non-perturbative gluon coming from the Vg q 
pairing above described can be estimated by using the recombination model. Actually, 
as the origin of the non-perturbative gluon is the recombination of a valon with an 
antiquark of the same flavor, the momentum distribution of intrinsic gluons in the 
\po g) fluctuation is simply given by 

g'^^'ix^Ql) = \a,\' Pp,{x,Ql) 

= l^iP dyyq{y) {y - x)vNq{y - x) , (9) 

where q{x) is given by eq. (0) and VNq{x) is the distribution of the g-flavored valon 
in the nucleon given by P| 

105 

VNq{x) = ^V^(l - ■ (10) 

Coefficients N, a, and (3; coming from the definition of the initial perturbative q, the 
recombination function and the valon-antiquark distribution in eqs. (0), and 
respectively; were included in the definition of |aip, the probability of the fluctuation. 
The intrinsic gluon probability distribution, Ppg{x,Ql), is accordingly normalized to 
unity. 

The intrinsic gluon probability density in the |A+ g) fluctuation can be estimated 
on similar basis. However, this fluctuation should be suppressed with respect to the 
\po g) one. 

Another way to have a proton fluctuation containing a proton-like object and 
an unflavored neutral meson would be through the self-recombination of the q — 
q pair produced by the gluon splitting of eq. (see Fig. In this case, the 
unflavored neutral meson must be a vector meson like a p*^ or u;. This is necessary to 
preserve the vectorial character of the initial perturbative gluon. However, this kind 
of fluctuations, consisting of a disconnected p° or uo and strongly suppressed 

by the OZI rule. 



6 



In Fig. ^, the intrinsic gluon distribution at the valon scale given by eq. is 
compared to the initial GRV-94 HO gluon distribution and the valence gluon 
distribution calculated in a Monte Carlo based model of the proton |Q 



3 Non-perturbative structure of pions 

Similar mechanisms should be at work in other physical hadrons, like pions. Indeed, 
if we expand the pion wave function as 



TT 



±,0 



bo TT^'O) + h fff'^g) + J2bmM:) , (11) 



we can identify the would be Ivr^'^vr^) fluctuation with the {Tft'^g) one, as we made for 
nucleons. In this way, the intrinsic gluon contribution to the pion low structure 
is given by a similar expression to that of eq. (^, 

(\ 2;) 

= l^iP / dyyq^iy) {y ~ x)vu{y - x) , (12) 

X Jo 

where qn{x) is an antiquark distribution analogous to q{x) in eq. (H) but with a valon 
distribution for pions vu{x) = 1, as given in [Q. Regarding the exponent related to 
velocity correlation, it turns out to be simply a = 1. For the valon scale in pions we 
have used the same value than for nucleons, Qf, ~ 0.64 GeV^. 

As explained above, a. Vqq pair should recombine into a gluon instead of forming 
a neutral pion structure. Thus the same conclusions drawn for nucleons hold for 
fluctuations of the pion into a Fock state containing the initial pion plus a neutral 
pion-like object. 

In Fig. ^, the g^^{x, Ql) predicted by the model is shown, and compared to the 



GRV-P HO distribution and the gluon distribution in pions obtained in Ref. 0. 
It is interesting to note that our model predicts intrinsic gluon distributions carrying 
more average momentum in pions than in nucleons. This is a consequence of the fact 
that the Vg q pair giving rise to the intrinsic gluon carries more average momentum 
in a I'n't'^ g) fluctuation than in a \po g) fluctuation. 

The fact that pions do not fluctuate into |7r^'°7r°) states but into \T^t'^ g) ones has 
additional consequences on their low scale structure. For example, for charged 

^In 1^, a model for hadrons is proposed in which primordial pdf corresponding to valence quarks 
and gluons are assumed to have Gaussian distributions with widths fixed from experimental data. 
These initial pdf are then complemented with contributions coming from \MB) fluctuations. In this 
model, intrinsic gluons are supposed to be present from the very beginning. In our model we are 
proposing a dynamical mechanism for their generation. This is the main difference between these 
two approaches. 
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pions the first contribution to the sum in the RHS of eq. (|TT]) arises from the 
{\K~K^)) fluctuation of the |7r+) (|7r~)) state. Thus, 
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(13) 
(14) 



Then the first contribution to the intrinsic qq sea arises in the strange sector and 
there are no uu and dd intrinsic seas in charged pions. The structure of charged pions 
at the low Ql scale is thus given by 



Vg/^± (X, Q 



&oP^^n(a;) + l&ip / —PgAy)'"u(- 

J-^ y Kvy 



+ \b 



SAX 



X y 



-PKKiy) VKq 



X 



Ql) = M' r^PKKiy)vKs (- 

Jx V \V. 



(15) 



where Vg/T,± are the resulting pion valence quark densities, vkq and vks are the light 
and strange valon distributions in Kaons, and Sg = Sg the non-perturbative strange 
quark distributions in charged pions. Pkk is the probability density of a Kaon inside a 
pion and Pg-^ix) = P7rg(l— x) is the charged pion distribution in the \Tf'^ g) fluctuation. 
The hadronic distributions inside pions, P^^g and Prk, are given by similar formulas 
to those of eqs. (0)-(|^). 

It is interesting to note that light quarks in a \KK) fluctuation contribute to 
the charged pion low valence densities but not to their intrinsic sea distributions. 
This is because although there are non-perturbative contributions to the light quark 
distributions in charged pions, they appear in the u (u) and d (d) sectors but not 
in the u{u) and d{d) for the vr" (vr^) respectively f\. In turn, for neutral pions the 
hadronic Fock state expansion has the form 



TT 



+ 61 7r°(7) + 6; 



+ 



K- K+) - 



V2 

Then, by analogy with eqs. (^), we can define 

hbo\'vnix) + l\b,\' f'^PgMvJ- 
2 2 Jx y \y^ 

'dy 



(16) 



Vg/nO{x,Ql] 



\b3 



X 



PKK{y)vKq , 

^ y \y 



•^Recall the flavor structure of tlie particles involved: 7r+(w(i) 
K-{us)K°{sd). 



(17) 

K+{us)K°{sd) and 7r-(t2d) ^ 
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for the valence quark densities at the low Ql scale, and 

Su/no {x, Ql) = (x, Ql) = (x, Ql) = (x, Ql) 

= \h\' f'^PUy)vn(-] 
y \yj 

Ss{x,Ql) = Ss{x,Ql) = Ihl'^ —PKK{y)vKs{-], (1 

■fo^ y \y 



for the intrinsic up, down and strange seas. Gluon distributions are given by eq. (|T2D 
for both neutral and charged pions. 

It should be noted that, although considering the Sg/^o densities as part of the 
intrinsic light sea or part of the valence densities in the 7r° is a matter of convention, 
the low Ql structure of the 7r° is different to the structure of charged pions. The 
difference is precisely given by the contribution of the |7r^7r+) fluctuation which can 
only occur in a tt" state. 

As a final result, notice that the intrinsic quark-antiquark sea of pions turns out 
to be symmetric as a consequence of the hadronic structure of the fluctuations. This 
is in contrast to the tipically unequal intrinsic quark and antiquark distributions of 
the nucleon (see e.g. [0] and Refs. therein). 



4 Conclusions 

In this paper we have analysed some important consequences of making a hadronic 
Fock state expansion of the nucleon and pion low wave-functions out of a novel 
mechanism for generating the cloud. We have shown that within such scheme it is 
possible to generate not only non-perturbative quark-antiquark distributions but also 
the gluon sea needed at the low starting scale for DGLAP evolution. 

These non-pertubative quarks and gluons are responsible for the bound nature 
of any hadron state, as they bring about the interactions between valence quarks. 
The non-perturbative quarks and gluons can be consistently identified with the so 
called intrinsic sea, in contrast to the extrinsic sea. On the other hand, the extrinsic 
sea is perturbatively generated by the probe momentum Q^, and is part of the own 
structure of valons, as discussed long time ago by Hwa [^]. 

In this sense our approach leads to a unification of two different pictures of the 
hadron structure; namely, the early picture of (non-interacting) valons f^, and the 
intrinsic sea idea of Brodsky et al. 0, which provides the binding agent for the 
bound hadron state. On the other hand, our approach allows a full representation 
of the non-perturbative processes giving rise to hadronic quantum fluctuations. This 
fluctuations are due to the perturbative production of a qq pair which recombines with 
the remaining valons. Thus a connection between the physics of hadronic reactions 
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and that of hadronic fluctuations is established through the well known recombination 
mechanism. 

A remarkable feature of the approach is that neutral pion fluctuations are here 
inhibited and, in turn, non-perturbative gluons take place. The reason is that neutral 
unflavored structures like the initial Vqq objects, are more likely to recombine rapidly 
into gluons than into neutral pions, in contrast to flavored structures like Vqq' which 
cannot do so. Thus, the hypothetical cloud of quantum fluctuations like |po 7r°) does 
not contribute to the sum over \BiMi) in the RHS of eq. (||) but to the second term, 
|po g), providing the source of valence-like gluons in the proton. 

Thus, within our scheme not only intrinsic quarks and antiquarks but also gluons 
are generated through quantum fluctuations of the low hadron ground state. 

Concerning pions, we have calculated their quark- ant iquark and gluon distribu- 
tions at low Q^. We have also shown that the non-perturbative structure of charged 
and neutral pions are different. The difference arises from the \tt^tt~) fluctuation 
appearing in the hadronic Fock state expansion of the 7r° wave-function but not in 
the charged pion ones. Finally, we have shown that the pionic intrinsic quark and 
antiquark distributions are symmetric, as a result of the specific features of its quan- 
tum fluctuations. This is in contrast to the structure of generic baryons which have 
asymmetric intrinsic quark and antiquark distributions 0, |T0[. However, it should be 
noted that for mesons containing a light and a heavier valence quark, the situation is 
different and the intrinsic quark- ant iquark sea must be asymmetric [|I^ . 

Summarizing, we have proposed a possible scenario for the origin of the valence- 
like sea quark and gluon distributions nedeed at the low (input) scale in order to 
describe the experimental DIS data for nucleons and pions. We have also discussed 
the low scale structure of charged pions and shown that, besides valence quarks, the 
model predicts only gluon and strange intrinsic sea distributions as a suitable low 
starting point for perturbative DGLAP evolution. On the other hand, for neutral 
pions, intrinsic light quark-antiquark distribution have to be considered as well. This 
signals a remarkable difference between the non-perturbative structure of neutral and 
charged pions. 
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Figure 1: Diagramatic representation of the process giving rise to a |7r~A"'"+) fluc- 
tuation of tlie proton (upper). Process leading to the generation of intrinsic gluons 
through the recombination oi avq q pair into a gluon (lower) . 
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Figure 3: Intrinsic gluon distribution for nucleons predicted by the model at Ql ~ 0.64 
GeV^ (solid line) compared to the valence gluon distribution given by the model of 
Ref. § at ~ 1 GeV^ (dashed line) and the initial GRV-94 HO gluon distribution 
at Ql = 0.4 GeV^ (point line). Model curves were normalized to the value of the 
integral over x of the GRV-HO gluon distribution. 
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Figure 4: Intrinsic gluon distribution for pions predicted by the model at Ql ~ 0.64 
GeV^ (solid line) compared to the valence gluon distribution given by the model of 
Ref. at some hadronic low Qq scale (dashed line) and the initial GRV-P HO gluon 
^0 = 0.3 GeV^ (point line). Model curves were normalized to the 



distribution 13 



at 



value of the integral over x of the GRV-P HO gluon distribution. 
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